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We calculate non-linear Knight Shift K vs. susceptibility x anomalies for Ce ions possessing local moments in metals. The 
ions are modeled with the Anderson Hamiltonian and studied within the non-crossing approximation (NCA). The K — vs. — x 
non-linearity diminishes with decreasing Kondo temperature To and nuclear spin- local moment separation. Treating the Ce 
ions as an incoherent array in CeSna, we find excellent agreement with the observed Sn K(T) data. 
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The origin of Knight shift anomalies in metals with 
localized moments that undergo the Kondo effect has 
been a subject of great interest in the condensed mat- 
ter community over a period of nearly twenty five years 
p]-pT|. The central physical concept is that the many 
body screening cloud surrounding a Kondo impurity site 
should give rise to an anomalous temperature dependent 
Knight shift at nuclear sites due to the coupling of the 
local moment to the nuclear spin through the screening 
cloud ll|J^-|5|,p|,p^ . This would be manifest in a devia- 
tion from a linear relation of the Knight shift K to the 
magnetic susceptibility x below the Kondo temperature 
Tq. Another way to say this is that in the absence of 
an anomaly, the contribution K^r, T) from a local mo- 
ment at distance r from the nucleus can be written as 
f{r)x{T). This factorization does not hold if there is 
an anomaly (instead K{r,T) = f{r,T)x{T) due to the 
temperature dependent polarization cloud). The clas- 
sic experiments by Boyce and Slichter on the low 
Kondo temperature {Tq ~ lOK) alloy Fe:Cu displayed 
no evidence for this anomalous Knight shift behavior. In 
contrast, pronounced Knight shift anomalies have been 
observed in the concentrated heavy electron materials 
CeSna [||| and YbCuAl |[||, which have been described 
as Kondo lattice systems with Tq ~ AOOK. In view of 
the Boyce- Slichter result, the question is raised whether 
these anomalies represent a coherent effect of the periodic 
lattice rather than a single ion effect. However, recent 
experiments on the proposed quadrupolar Kondo alloy 
[ p2| , p^ Yi-a-U^Pda demonstrate that for concentrations 
of 0.1-0.2 there are pronounced non-linearities in the Y 
Knight shift for sufficiently large distances away from the 
U ion fill. 

In this paper, we present some results from our sys- 
tematic theoretical studies of the Knight shift in heavy 
electron materials. We find, as expected from earlier an- 
alytic theories, that the magnitude of the Knight shift 
anomaly is reduced with decreasing Kondo temperature 
and decreasing distance between local moment and impu- 



rity sites. Ours is the first study to display this effect in 
a realistic model calculation. We find that the Knight 
shift anomaly of CeSna can be accounted for by inco- 
herent, single ion physics. In particular, by carrying out 
a full lattice sum on CeSns, we obtain an excellent one 
parameter fit of the calculated Knight shift to the ex- 
perimental one, with a small low temperature deviation 
possibly due to the over-simplifications of our model. We 
conjecture that coherence effects may explain some of the 
small numerical disagreement. 

The application of NMR in Kondo systems received 
impetus from Heeger ^ who suggested that the anoma- 
lous cloud was detectable in Knight shift measurements. 
Essentially, the oscillatory conduction spin density s{r) 
induces a local moment-nuclear moment interaction. In 
the limit Tq — > 0, this interaction X(r) is well described by 
the second order perturbation theory Ruderman-Kittel- 
Kasuya-Yosida expression (s(r) ^ T{r) ~ cos{2kFr)/r^). 
Ishii confirmed that for an S*/ = 1/2 pure-spin 

Kondo impurity (with no charge fluctuations) an anoma- 
lous conduction spin density cloud sets in beyond the 
Kondo screening length S^k — ^vp/ksTo, where vp is 
the Fermi velocity. Inside this radius, conventional tem- 
perature independent RKKY oscillations dominate of the 
kind observed by Boyce and Slichter. Outside the screen- 
ing length, at T = 0, the anomalous term will dominate 
also with an RKKY form but an amplitude of order D, 
the conduction bandwidth, compared with D{N{0)J)'^ 
for the Ruderman-Kittel term, where A'^(O) is the con- 
duction electron density of states at the Fermi energy 
and J the conduction electron- local moment exchange 
coupling. Ishii did not calculate the explicit tempera- 
ture dependence of this structure, but did anticipate that 
it would vanish above the Kondo scale. Scaling analy- 
sis confirmed the asymptotic factorization of the Knight 
shift for short distance and low temperature j|] . A possi- 
ble understanding of the Boyce and Slichter results, then, 
is the Cu nuclei they sampled were at distances r << S.K 
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from the Fe ions. 

In an Anderson model treatment of the problem, 
charge fluctuations are allowed. Consider, for example, 
a Ce ion with nominally one / electron giving rise to 
the local moment. The ensuing Kondo effect is best un- 
derstood in terms of a mixing of the empty configura- 
tion with one in which the Ce moment is screened by a 
superposition of conduction hole states. This opens up 
new possibilities for detecting the anomalous screening 
cloud. One interpretation of the Boyce-Slichter results is 
that due to subtleties of spin conservation in the singlet 
ground state the spin cloud is unobservable, while the 
anomalous charge cloud is observable through electronic 
field gradient measurements at atomic nuclei ||]. 

We shall focus on the second possibility, that the 
anomalous spin cloud now has a finite amplitude within 
r < £,K which is proportional to the weight of /° in the 
ground state. This anomaly scales approximately lin- 
early with Tq, so the cloud will be virtually undetectable 
in small Tq systems, but may be observable in high Tq 
systems. Our results support this picture as we shall 
explain in more detail. 

Our starting point is an Anderson model Hamiltonian. 
We shall first discuss the situation for Ce^+ and Yb'^^ 
ions, and write down the model only for the Ce case. For 
a single Ce site at the origin, the model is 

n = Hc+Hf+Hcf+H^ (f) 

with 

ka 

the conduction band term for electrons with a broad 
featureless density of states of width D, taken to be 
Lorentzian here for convenience; with 

nf = J2^n\fj^><fM (3) 

jm 

where j — 5/2, 7/2 indexes the angular momentum mul- 
tiplets of the Ce ion having azimuthal quantum numbers 
m, with e/5/2 = -2eV,efT/2 = e/5/2 + Aso = -l.7eV (we 
take the /" configuration at zero energy); with 

kjma 

where Vj:^^^ = VY^^.^ik) <3m-a, l/2a\jm > ///V7, 
V being the one particle hybridization strength and Ng 
the number of sites; with 

Hz = -HBH42^anj: ^-^gjm\fjm >< fjm\] 

kcr 

(5) 



being the Zeeman energy of the electronic system for a 
magnetic field Hz applied along the z-axis. In addition to 
this, we must add a term coupling the nuclear spin system 
to the conduction electrons, which we take to be of a 
simple contact form ~ vecI{f)-S(fi) for each nuclear spin 
I{r) at position r with S{R) the conduction spin density 
at the nuclear site, and a nuclear Zeeman term. In terms 
of the parameters, the Kondo scale characterizing the low 
energy physics is given by 

kBTQ = D{^f'\^^f'^eM"-^) (6) 

where the hybridization width F = tiN{Q)V'^ . 

We treat the Anderson hamiltonian with the non- 
crossing approximation (NCA), a self- consistent dia- 
grammatic perturbation theory discussed at length in the 
paper of Bickers et al. Jl4[ . This approximation is con- 
trolled by the large orbital degeneracy of the Ce ground 
state. It does show pathological behavior (due to the 
truncation of the diagrammatic expansion) for a tem- 
perature scale Tp « Tq in this conventional Anderson 
model, provided the occupancy > 0.7. In prac- 
tice, this is not a problem for iV > 4 as shown in Ref. 
[Q , in that comparison of NCA results with exact ther- 
modynamics from the Bethe-Ansatz shows agreement at 
the few percent level above Tp. Hence, this is a reliable 
method for our purposes. 

The approach in the NCA is to write a propagator 
for each ionic state of the Ce site (i.e., f^j = 5/2,7/2 
and /° in the present model), solve self-consistent inte- 
gral equations to second order in the hybridization for 
the ionic propagator self-energies, and then to calculate 
physical properties as convolutions of these propagators. 
To evaluate the Knight shift we employ the lowest order 
diagram coupling nuclear spin to Ce magnetic moment, 
as shown in Fig. I. The full expression corresponding 
to this diagram is cumbersome and shall be presented 
in detail elsewhere. To the extent that the dynamics 
of the empty orbital can be neglected, this expression 
factorizes into a nearly temperature independent RKKY 
interaction (modified due to the spin-orbit coupling and 
anisotropic hybridization from the original form) times 
the /-electron susceptibility. Thus, no anomaly results 
from the diagram in this limit. In this limit, the sus- 
ceptibility in the diagram corresponds to the the leading 
order estimate used in Ref. to compare with exact 
Bethe-Ansatz results. 

However, for T — 0, the empty orbital propagator may 
be written in an approximate two-pole form, one with 
amplitude 1 — Z , Z = Trfc^To/GF, centered near zero 
energy, and one with amplitude Z centered at ~ 6^5/2 — 
ksTo which reflects the anomalous ground state mixing 
due to the Kondo effect. The first term gives conventional 
RKKY oscillations modulo the anisotropy and altered 
range dependence induced by the m, k dependence of the 
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hybridization. The amphtude of the second term goes to 
zero above the Kondo temperature. It is this term which 
may be traced to the anticipated anomalous Knight shift, 
and within such a two-pole approximation may be seen to 
be finite within £^k, have a stronger distance dependence 
in that regime, but possess an amplitude of order Z only 
within this distance regime. Beyond the amplitude 
is of order 1 /N and the shape of the spin oscillations is 
the same as that found from the high frequency pole of 
the empty orbital propagator. 

The diagram of Fig. 1 has been studied with the NCA 
previously pO[|, but only for the spin 1/2 model with 
infinite Coulomb repulsion, and for a limited parameter 
regime (very low Tq values) and short distances (r << 
£,k)- In consequence, no strong evidence was found for a 
Knight shift anomaly in this previous work. 

Our numerical procedure, briefly, consists of solving 
the NCA integral for the Anderson Hamiltonian speci- 
fied above on a logarithmic mesh with order 600 points 
chosen to be centered about the singular structures near 
the ground state energy Eo ^ e/5/2- We then feed the 
self-consistent propagators for the empty and doubly oc- 
cupied orbitals into the convolution integrals obtained 
from the diagram of Fig. 1, which allows for evalua- 
tion of the Knight shift at arbitrary angle and distance 
from the nuclear site. It is convenient to take the nuclear 
site as the origin in this case leading to phase factors 
g-ifc _R -j-j^g hybridization Hamiltonian Ti-cf, where r is 
the nuclear-Ce site separation. These factors give the 
oscillations and position space angular dependence in K. 
All contributions from the Ce susceptibility are included. 

To demonstrate the dependence of the K — vs. — x 
anomaly as a function of Kondo scale Tq, we have eval- 
uated the diagram of Fig. 1 for a single local moment 
placed at a fixed distance r = 3.3kp^ — 3.3A for our 
choice of kp — lA^^ and angle of 0° with respect the 
quantization axis. We have tuned the Kondo scale by 
holding fixed all parameters except the hybridization. 
The Knight shift is scaled to a susceptibility by match- 
ing at high temperatures, and the susceptibility units are 
scaled by a fit of one calculation to the data for CeSna 
assuming D — 3eV. The results are shown in Fig. 2(a). 
Clearly, as the Kondo scale is reduced, the magnitude of 
the deviation between a linear Kvs.x relation is system- 
atically reduced. Also, as the separation r is decreased 
(c.f., Fig. 2(b)) the magnitude of the K — vs. — x non- 
linearity diminishes. Very similar results are obtained 
for Yb compounds (the Yb'^"'" ion has a lone 4f hole and 
our procedure decribes these with a simple particle hole 
transformation, which we shall discuss in detail in a sub- 
sequent publication). 

To assess the relevance of this single site physics to the 
periodic compound CeSna we have performed a lattice 
sum about a given Sn nucleus from all surrounding Ce 
ion. We assume the Knight shift contribution of each ion 



to be described by this single site physics, known to be 
a good approximation at high temperatures where the 
ions are incoherent with one another, and known to pro- 
vide a very accurate description of the thermodynamics 
in many cases. We carried out the sum to several hun- 
dred shells about the Sn site, obtaining good convergence 
at all calculated temperatures. We fixed the parameters 
by fitting the susceptibility data of CeSns. The result 
is shown in Fig. 3, where we have scaled the Knight 
shift by an intermediate range temperature match to the 
susceptibility (note that the NMR data of Ref. |^ ex- 
tends only to room temperature). There are several no- 
table features here: (i) the calculated amplitude of our 
Knight shift prior to scaling is actually negative, which 
implies that the fit is sensible only if the assumed contact 
coupling between conduction and nuclear spins is nega- 
tive, which actually makes sense because the Sn nucleus 
should dominantly couple through core polarization; (ii) 
the magnitude of the anomaly goes in the right direction 
and begins at the right temperature to agree with the ex- 
perimental anomaly , though the quantitative value is 
slightly too high-given the highly oversimplified conduc- 
tion band we are employing, we don't regard this as a se- 
rious defect in the calculation; (iii) the magnitude of the 
anomalous contribution from the Ce sites does go down 
with distance from the Sn nucleus-the theoretical data 
which most closely match those of experiment are taken 
at the fixed distance r = 3k = 3 A and angle of 0°. 
However, the anomaly is then still surprisingly large even 
at this short distance. While we cannot rule out that the 
discrepancy reflects our oversimplified band structure, we 
conjecture that lattice coherence effects may drive our 
calculated lattice sum in this direction, because the re- 
peated scattering of conduction electrons off of the Ce 
ions (now in every unit cell) has the effect of reducing 
the screening length We intend to study this idea 

further within a local approximation {d = oo expansion) 
to the lattice model. We note that the high temperature 
deviation of Knight shift and susceptibility may be real- 
istic and is related to the different couplings to ground 
and excited spin orbit multiplets. 

In summary, we have computed Knight shift anomalies 
within a realistic model for Ce ions in metals for the first 
time. We find that the magnitude of the non-linearity 
scales down in size as the Kondo scale is diminished or 
the nuclear moment-local moment separation is reduced. 
Modeling the Ce ions incoherently in CeSna, we find that 
summing over all nattice sites the resulting Sn Knight 
shift agrees quite well with experiment. The remaining 
quantitative discrepancies may be due to a combination 
of our oversimplified band structure for the conduction 
states together with lattice coherence effects. 
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Figure 1. Feynnian diagram of coupling between Ce 
local moment and nuclear spin in the infinite U Anderson 
model. The dashed line represents the singly occupied 
(/^) state and wavy line the empty orbital (/°) state. 
The solid line represents the conduction electron. The 
dot-dashed line is a propagator for the nuclear spin states 
and H is the extenal magnetic field. 

Figure 2 (a) Calculated Knight shift K{T) vs. suscep- 
tibility x(T) for a single Ce site at kpT = 3.33 from a nu- 
clear moment and angle = 0. Fixing the /-level energy 
65/2 = — 2eV, and the spin orbit splitting Ago = Q.ZeV, 
the hybridization is varied to illustrate the dependence of 
the nonlinearity on the magnitude of Tq. The diagram of 
Fig. 1 is used to calculate K{T). The magnitude of the 
non-linearity diminishes as Tq is reduced. The theoreti- 
cal Knight shift has been shifted by a common offset and 
scale factor to match the susceptibility, (b) Calculated 
Knight shift K{T) for a single Ce site vs. susceptiblity 
X(T) for varying separation with the Kondo scale used to 



fit the CeSna x(T) (see Fig. 3). For each plot the angle 
is held at 9 with respect to the nuclear moment-Ce axis. 
The magnitude of the non-linearity diminishes as kpT is 
reduced. The theoretical Knight shifts have been shifted 
by offset and scale factors to match the susceptibilities; 
this does not affect the relative size of the anomaly. 

Figure 3 Temperature dependence of Sn Knight Shift 
K{T) and Ce x{T) (both calculated and experimental 
results Q) for CeSna. The theoretical K{T) is calculated 
using the diagram of Fig. 1 and with Tq chosen to fit the 
experimental x{T) data. A full (incoherent) lattice sum 
is carried out over several hundred shells of atoms. 
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